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Abstract 

Q ' In the present paper, we study the scalar sector of the quasinormal modes of 

^D . charged general relativistic, static, and spherically symmetric black holes coupled 

to nonlinear electrodynamics and embedded in a class of scalar-tensor theories. 
We find that for certain domain of the parametric space, there exist unstable 
^ ■ quasinormal modes. The presence of instabilities implies the existence of scalar- 

H I tensor black holes with primary hair that bifurcate from the embedded general 

relativistic black-hole solutions at critical values of the parameters corresponding 
to the static zero-modes. We prove that such scalar-tensor black holes really 
exist by solving the full system of scalar-tensor field equations for the static, 
spherically symmetric case. The obtained solutions for the hairy black holes 
are non-unique, and they are in one to one correspondence with the bounded 
states of the potential governing the linear perturbations of the scalar field. The 
stability of the non-unique hairy black holes is also examined, and we find that 
the solutions for which the scalar field has zeros are unstable against radial 
perturbations. The paper ends with a discussion of possible formulations of a 
new classification conjecture. 



*E-mail: ddoneva@phys.uni-sofia.bg 
■fE-mail: yazad@phys.uni-sofla.bg 
■l-E-mail: kostas.kokkotas@uni-tuebingen.de 
''E-mail: izhivkov@yahoo.com 



PACS numbers: 04.70.Bw, 04.50.Kd, 04.25.D- 



1 Introduction 

The quasinormal modes of black holes were intensively studied during the last decade. 
The great interest in the quasinormal modes is motivated on the one side by astro- 
physics and especially by the expected recent observation of one of the most important 
predictions of general relativity - the gravitational waves. Information about the phys- 
ical properties of the black holes (and other compact objects) can be obtained through 
the emission of gravitational waves and the characteristic frequencies of ringing, namely 
the frequencies of the quasi- normal modes (QNMs) [H [2]. The QNMs would allow us 
to distinguish between different compact objects - black holes, stars, to distinguish 
between different theories of gravity since they predict different characteristic spectra, 
to determine the global asymptotic charges like mass, charge, and angular momentum 
of the observed black holes and so on ^ . 

On the other side, the quasinormal modes are interesting and important for intrinsic 
theoretical reasons. During the last decade, the quasinormal modes became a powerful 
tool for the study of various pure theoretical problems which deepen our knowledge of 
spacetime and gravity, especially in the regime of strong fields j4j. They may serve, for 
example, as an indicator of black-hole phase transitions [5]-|10j. 

In the present paper, we study the scalar sector of the quasinormal modes of general 
relativistic, static, and spherically symmetric black-hole solutions coupled to nonlin- 
ear electrodynamics which are also solutions to the field equations of a certain class 
of scalar-tensor theories with nonlinear electrodynamics, i.e., they can be viewed as 
embedded in the scalar-tensor theories under consideration. In this paper, we consider 
a class of scalar-tensor theories defined by zero potential V{ip) = and a coupling 
function 

A{ip) = e-2^^' (1) 

corresponding to a{ip) = /Sip where /3 is a constanl|l|. These scalar-tensor theoriecl are 
indistinguishable from the general relativity in weak field regime (the post-Newtonian 
parameters have the same values as in GR), but can differ seriously from it in the strong 
field regime as it was first shown in [11] for neutron stars. It is natural to try to find 
similar non-perturbative deviation from GR when black holes are considered as sources 
of strong gravitational filed. According to several no- scalar-hair theorems, however, in 
the case when no matter sources are present in the theory [12] and in the case when the 
scalar-tensor theories are coupled to linear electrodynamics [13], the static, spherically 
symmetric black-hole solutions are indistinguishable from the corresponding solutions 
in GR (see also [H] - [IZ]). These theorems can be eluded when the scalar-tensor 
theories are coupled to non-linear electrodynamics since it has a non-vanishing trace of 



^The general definition of the function a{ip) is a{ip) = dlnA(ip)/dip and more details can be found 
in the next section. 

^ In the present paper we fix the cosmological value of the scalar field to be ipoo = 0. 



the energy-momentum tensor, and this nonzero trace acts as a source of the scalar field 
giving rise to black holes with nontrivial scalar field. Such scalar-tensor black holes 
were first found in [18j and [19j. Within the particular class of scalar-tensor theories 
considered here, the static and spherically symmetric black holes coupled to nonlinear 
Born-Infeld electrodjTiamics were studied in [20|. It was shown that black holes with 
nontrivial scalar field can exist only for /3 < 0. The most interesting result in [20] is the 
fact that the considered class of scalar tensor theories with /3 < allows the existence 
of non-unique black-hole solutions with the same mass and charge. This disproves the 
no-hair conjecture according to which the static scalar-tensor black holes are uniquely 
specified by their conserved asymptotic charges, namely, the mass and the charge. In 
other words, the scalar-tensor black holes can support primary scalar hair in contrast 
with the secondary scalar hair uniquely determined by the mass and the charge. 

Among the non-unique solutions present in the considered theory is the solution 
with zero (trivial) scalar field which is just the solution of the pure Einstein equations 
coupled to nonlinear Born-Infeld electrodynamics. From now on, the embedding of 
this general relativistic solution in the scalar-tensor theories under consideration will be 
called "t/ie trivial solution" . It turns out that in a certain class of scalar-tensor theories, 
the perturbations of the trivial solution form two independent sectors - scalar and 
metric- vector sector. In other words, the perturbations of the scalar field decouple from 
the perturbations of the metric and electromagnetic field. The equations governing 
the perturbations of the metric and the electromagnetic field coincide with those of 
pure general relativity. Since we are interested in the scalar-tensor theories, we focus 
only on the scalar sector. The equation governing the scalar perturbations feels not 
only the background metric and background electromagnetic field but also involves a 
term coming from the scalar-tensor theory under consideration. The presence of this 
scalar-tensor term in the potential leads to interesting consequences. Especially, the 
scalar-tensor term with /3 < gives negative contribution to the potential of the scalar 
perturbations and a negative minimum in the potential appear for a certain domain 
of the parameter space. The bound states associated with the negative minimum of 
the potential correspond to unstable modes. Therefore, for a certain sector of the 
parameter space, the trivial solution becomes unstable. How the instability of the 
trivial solution should be interpreted? Our interpretation is that the instability is a 
sign for the existence of new black-hole solutions with nontrivial primary scalar hair 
which bifurcate from the trivial solution at critical parameters corresponding to the 
static zero-modes. We prove this interpretation by numerically solving the full system 
of static, spherically symmetric scalar-tensor field equations and finding the explicit 
solutions with nontrivial scalar field. The number of nontrivial scalar-tensor black-hole 
solutions that bifurcate from the trivial one turns out to be in one to one correspondence 
with the number of bounded states of the scalar perturbations potential. In this way, 
we find a wide spectrum of new solutions which is much richer than what was previously 
found in [20]. 

We examine also the stability of the nontrivial solutions. We show that the non- 
trivial solutions for which the scalar field has zeros are unstable against radial pertur- 
bations. 

The paper is organized as follows. In section 2, we introduce the necessary back- 



ground of general equations. In section 3, the equation governing the scalar pertur- 
bations is derived. The numerical methods for the calculations of the QNMs and 
the numerical results are presented in section 4. The spherically symmetric solutions 
describing hairy black holes and bifurcating from the trivial solution are found numer- 
ically in section 5. In section 6, we examine the stability of the nontrivial black-hole 
solutions. The paper ends with discussion. 

2 General equations 

The action of the scalar-tensor theories coupled to nonlinear electrodynamics in the 
Einstein frame is 



-^J d'xV^A\v)L{X,Y), (2) 



where G* is the bare gravitational constant, g is the determinant of the Einstein frame 
metric (7^,^, R is the Ricci scalar curvature with respect to the metric (7^1^, ip is the scalar 
field, V{ip) is the potential of the scalar field, A{ip) is the coupling function between the 
sources of gravity and the scalar field and L{X, Y) is the Lagrangian of the nonlinear 
electrodynamics. The functions X and Y are given by 

X = ^^F,^F^^, (3) 

Y = :^^F,,{^Fr, (4) 

where F^^, is the electromagnetic field strength tensor and "V stands for the Hodge 
dual with respect to the metric g^^. 

The field equations derived from the action are given by 

R,, = 2d,vd,v + 2V{v)g,, - 2dxL{X, Y) (f^^F^^ - ^g.^F^t^F'^A 
-2A\^) [L{X, Y) - YdyLiX, Y)] g,,, 



V, [dxLiX, Y)F^-' + dyLiX, r)(^F)n = 0, (5) 

V^V^V^ = ^^ - MvH\^) [HX, Y) - XdxHX, Y) - YdyLiX, Y)] , 



where 

dlnA{^) , . 

a{^) = -^^. (6) 



In the present paper, we consider a class of scalar-tensor theories defined by zero 
potential V{ip) = and a coupling function 



A^) 
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/3^2 



(7) 



corresponding to a{ip) = Pip where /9 is a constant. According to some estimations the 
parameter /3 must obey the inequality P > —4.5 in order for the theory to be consistent 
with the experimental data [21]. That is why we shall restrict our considerations in 
this region of values of (3. The scalar-tensor theories with a{(p) = (3ip (and cosmological 
value (foQ = 0) as we have already mentioned are indistinguishable from general rela- 
tivity in the weak field regime but can exhibit non-perturbative effects and can differ 
seriously from general relativity in the strong field regime [llj.^20j. 

In order to be specific, we will consider here the Born-Infeld nonlinear electrody- 
namics. An important property of this nonlinear electrodynamics is that it is invariant 
under electric-magnetic duality rotation so it is enough to consider only magnetically 
charged solutions, and the electrically charged can be simply obtained from them. For 
pure magnetic solutions, Y = and the truncated Born-Infeld Lagrangian is 



L{X) = 26 1 




(8) 



where 6 is a parameter and in the limit 6 — )■ oo the Maxwell electrodynamics is restored. 
It can be easily shown that for the Born-Infeld theory, the following inequality holds 



XdxL{X) - L{X) > 0. 



(9) 



We will consider static, spherically symmetric, asymptotically flat black holes and 
the ansatz for the metric is 

ds"" = g.^dx^dx" = -f(r)e-^^^'Ut^ + ttt + ^' (dO^ + sin' ed(f)^) . (10) 

f{r) 

In this case, the function X reduces to 



X 



A-\ip) P2 



„4 ' 



(11) 



where P is the magnetic charge [20l[T8]. Respectively, the fleld equations are reduced 
to the following system of ordinary differential equations 



dS f dip 

dr V dr 



dm 

dr 

d_ 
dr 



i;(^) _^,„.„,) 



dr J 



\dr J 



(12) 
(13) 



- V7 ) = -M -4«(^)^'(V') [L - XdxL{X)] -rf(^\ ), (14) 



where we have introduced the new function ■m(r) defined by /(r) = 1 — 2m{r)/r. 
Without loss of generahty we consider scalar fields that vanish at infinity. Taking into 
account that a{!f) = f3ip it is not difficult to see that every static, spherically symmetric 
solution of the Einstein equations coupled to a nonlinear electrodynamics is also a 
solution to the above system with ip{r) = 0. The general relativistic solution describing 
black hole coupled to a nonlinear electrodynamics is naturally embedded in the class 
of scalar-tensor theories under consideration and, as we have already mentioned, will 
be called the trivial solution. The metric functions of the trivial solution are given by 

6{r) = 



2M 4br^ / / ^2 , 

/W = i-- + — (i-f + a^Fil '^^> 

4P2 /I 1 5 p2 



3r2 V4'2'4' 26H 

where J-" is the hypergeometric function |22] and M is the black hole mass. 

For the trivial solution, it can be easily proven that extremal black-hole solutions 
can exist only for bP^ > 1/8. However, no extremal black holes can exists in the 
considered class of scalar-tensor theories when the scalar field is nontriviao 



3 Perturbations of the trivial solution 

In the present paper, we consider the perturbations of the trivial solution within the 
framework of the described class of scalar-tensor theories. The metric perturbations 
will be denoted by 5g^,y, the electromagnetic perturbation^l by ^A^ and the scalar ones 
by 6ip. Since the trivial solution has zero scalar field, it is not difficult to see that in 
the considered class of scalar-tensor theories the equations governing the perturbations 
of the metric 6g^u and the electromagnetic field SAfj^ are decoupled from the equation 
governing the perturbations 6(f of the scalar field. The equations for 6gf^,y and 6Af^ are 
in fact the same as those in the pure Einstein gravity. Since we are interested in black 
holes in the scalar-tensor theories we will focus only on the scalar perturbations S(p. 

From the third field equation in ([5]), we obtain the following equation governing the 
perturbations of the scalar field 

V(°)V(°)'^5^ = 4/3[X(°)a^,o,L (X(°)) - L (XW)]5^, (16) 

where the covariant derivative V/j and the function X^^^ are with respect to the un- 
perturbed background solution. Since the background solution is static and spherically 
symmetric, the variables can be separated in the following way 

5^ = ^e^^%U0,^), (17) 



^A detailed proof can be found in [50] 

^Here, A^ is the gauge potential, i.e., F^,^ = d^Ai, — d^Afj 



where a; is a parameter and Yim(0,4>) are the spherical harmonics. After substituting 
in fll6p we find 

d f j.dip\ ^ I- 2 



^Ud^j + l^-^**'^''- (^**) 



where the potential U is defined by 



f/(r) = / 



-f + ^^ + 4/3[xW9^(o)L (XW) - L (XW)]' 



The equation (1181) can be cast in the Schrodinger form 

d'^ilj{r 



(19) 



^, +uj'^{n) = U{n)^{n) (20) 

dr^ 
by introducing the tortoise coordinate 

which maps the domain r G [r^^, oo) {th is the radius of the horizon of the black hole) 

to r^ G (— oo, oo). 



4 Quasi-normal modes 

4.1 Dimensionless quantities 

The reduced system of equations (fT2l) - (fT^ is invariant under the rigid rescaling 

r -^ Ar, m -)■ Am, b -^ \'%, P -^ \P, (22) 

where A G (0, oo). Therefore, one may generate in this way a family of new solutions. 
The mass, the temperature, and the entropy of the new solutions are given by the 
formulas 

M^AM, T^X-^T, S^X'^S. (23) 



Respectively, the equation for the perturbations of the scalar field ( 1181) is invariant 
under the rigid rescaling 

r -)■ Ar, m -)■ Am, b -^ X^^b, P -)■ AP, u -)■ A^-^w. (24) 

As a consequence of the rigid symmetry we may restrict our study to the case 
b = 0.01 by introducing the dimensionless quantities 

T" 771 I fvl 

r^-, m-^— , 6->6e^ ^ ^ T' ^ ^ ^^. M ^ — , (25) 



where the choicqj 6 = 0.01 requires 



^^7^- 



'The very choice of the value 0.01 is for numerical convenience. 



4.2 Qualitative considerations and numerical algorithm 

Our task is to calculate the QNMs of the trivial solution embedded in the considered 
class of scalar-tensor theories. Thus we have to find the QNM frequencies governed by 
eq. flTSjl (or equivalently by eq. fl20|) ) where the boundary condition are the standard 
ones: purely outgoing wave at infinity and purely ingoing wave at the black hole hori- 
zon. We have to note that the problem reduces to studying the scalar perturbations of 
the pure Einstein-Born- Infeld black holes when the parameter in the coupling function 
/3 = and this case is studied in [23] |j. 

It is the potential U that determines the qualitative and quantitative behavior of the 
QNMs for the already fixed boundary conditions. That is why we will first comment 
on the qualitative behavior of U defined by eq. (fTOll . It can be easily checked that the 
only difference between the equation governing the perturbations of the scalar field (ITSll 
and the equation governing the scalar perturbations in the case of pure Einstein-Born- 
Infeld theory is the term Al3[X^^^dx(o)L (X(°)) - L (X^o))] in the potential ((19]). When 
we take into account equation ([9]), it is clear that the sign of this term is the same as 
the sign of /3. Thus, the potential can become negative when /3 < 0. On Fig. (H the 
potential U{r) for M = 1.0, /3 = —4 and / = is plotted for a range of parameters. As 
it can be seen, a positive maximum is present for all curves, and the numerical results 
show that such maximum also exists for all of the studied values of the parameters M, 
P, /3, and /. For large values of the charge, an additional negative minimum appears 
when /3 < 0, and it can lead to instabilities. When /3 > 0, the potential is qualitatively 
the same as for the pure Einstein case. 

For fixed values of the parameters, the depth of the potential well, if it exists, 
decreases with the increase of / because of the contribution of the second term in flTOl) . 
For high enough /, the potential becomes positive. 

The calculation of the QNMs in practice, however, faces serious numerical difficul- 
ties. The first reason for that is the complexity of the background solution which can 
not be expressed in elementary functions as well as the presence of negative minimum 
in the potential for some values of the parameters. This makes us use direct integration 
methods which has limited abilities since the more sophisticated and accurate methods 
are unapplicable. The second reason is the numerical instabilities that appear in the 
region of transition from stable to unstable modes which will be discussed below. It 
is well known that even for the Schwarzschild solution, the direct integration meth- 
ods face serious difficulties when the imaginary part of the quasinormal frequencies 
becomes greater than the real part. 

It turns out that the properties of the perturbation functions and the modes in 
the stable and unstable sector are different, and that is why it is useful to consider 
the stable and the unstable modes separately. More precisely, as we shall see below, 
the unstable modes define a self-adjoined Sturm-Liouville boundary value problem for 
which the eigenvalues are real. This simplifies the treatment of the problem and makes 
the numerical procedure stable. 



^Our results for /3 = are qualitatively similar to the results obtained by Fernando and Holbrook 
[23] but the computed numerical values of the frequencies are different. 
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M=1.0,/3^-4.0,/=0 




Figure 1: The potential f/ as a function of r for values of the parameters M = 1.0, 
/5 = —4.0, / = and for several values of the magnetic charge P. The range of values 
of P is chosen so that the formation of a negative minimum can be seen as we increase 
the charge. 



4.2.1 Stable modes 

The required QNMs boundary conditions lead to an exponential growth of the radial 
perturbation function for the stable modes as r^, — t- ±oo . To circumvent this difficulty, 
we use the numerical procedure introduced by Chandrasekhar and Detweiler [21], [25]. 
The idea of the method is to introduce a new function 0(r,,) defined by 

^ = 6'^^'^'^', (27) 

where ip is the radial part of the perturbations of the scalar field. Equation ( l20l) 
transforms into a Riccati type equation 






and the boundary conditions, according to the sign convention in (TT7|) . are 



— !■ CJ, 



-)■ -W. 



(28) 



(29) 



We apply the following shooting procedure for calculating the QNM frequencies. 
Equation (!28|) is integrated for test values of the complex parameter u and a QNM 
frequency u is found when the boundary conditions (12^ are satisfied within the re- 
quired accuracy. For each test value of uj two integrations are performed - inward and 
outward, and they are described below in details. 

First, we integrate equation f l28|) outward from a chosen point r=i,_oo to a matching 
point r^.match which is close to the maximum of the potential. The initial condition at 
r*_oo is computed using the expansion of the function %Ij at the black hole horizon 



'i/'i, 



Eft 



rH] 



(30) 



where r^ is the radius of the horizon and (3^ are constants (/3i = 1 and /S^ = for 
s < 1). The coefficients Ps {s > 1) can be calculated from the recurrence relation that 
is obtained when we substitute (15U]1 into eq. (12U1) . 

The second inward integration is from a chosen point r*oo III to the same match- 
ing point rematch and the initial condition at r*oo is calculated from the asymptotic 
expansion at infinity 

oo 

V'|,.^oo = e-*'^^* J]ay-^ (31) 

where as are constants (ai = 1 and a^ = for s < 1). The coefficients as (s > 1) can 
be found from a recurrence relation derived in the same way as for the Ps ■ 

After the two integrations, we compare the values of the function 0, resulting from 
the inward and outward integration, at the matching point r=nnatch, and a QNM fre- 
quency u is found when these values coincide within the required accuracy. 

As it is pointed out in [21], the numerical method we use may become unstable 
when the imaginary part of the QNM frequency becomes greater that the real one, and 
this situation actually occurs in our calculations in the region of transition from stable 
to unstable modes. That is why we have to use a stiff ODE solver, and during our 
tests, it turned out that the semi-implicit extrapolation method described in Numerical 
Recipes [26] is a very good choice. The search of QNM frequencies cu described above, 
can be considered as a root search procedure, and we used the Muller method for it 

4.2.2 Unstable modes 

According to our sign convention in flT7|) . the unstable QNM frequencies should have 
negative imaginary part. In this case, we can substitute u = ujr + iuj, where ur and 
ui are both real and uj is negative. Hence, the boundary conditions for the unstable 
quasi-normal modes are vanishing on both ends 

V'lr.^-oo ^ e^'^«'-*e-"^^* |,,^_oo ^ 0, (32) 

^Ir.^oo ^ e-'"«'-*e"^'-*|,,^oo -> 0. (33) 

But in this case the boundary value problem (BVP) discussed above is self-adjoint and 
has real eigenvalues u"^. Therefore, the unstable QNM frequencies should be purely 
imaginary, and they correspond to the bound states of the potential well. More infor- 
mation on this subject can be found in [27], and similar considerations, but for neutron 
starts in scalar-tensor theories, are made in [28] . 

The perturbation function ip for the unstable modes will have no zeros for the 
first bound state, one zero for the second bound state, and so on. The function cf) 
defined by fl27|) is divergent in the points where ip is equal to zero which makes the 
numerical procedure used by Chandrasekhar and Detweiler inapplicable for calculating 
the unstable modes (except for the first bound state where tp has no zeros). Hence, we 
have to solve the original equation for the ip function (120!) and the same expansions on 
the boundaries (l30!l and (13T|) can be used for calculating the boundary conditions. 



'''The point r^oo is chosen so that the ratio r^oo/li^l is kept fixed. 

10 



The numerical method used is quite similar to the shooting procedure described 
above: equation ( 120 p is integrated inward and outward for test values of the parameter 
a;, and because of the arbitrary linear scale, it is enough to match at the chosen point 
rematch only the quantities 

resulting from the two integrations. 

4.3 Results 

We calculate the QNM frequencies of the trivial solutions using the numerical tech- 
niques described above. As we said, we will restrict ourselves to calculating the QNM 
frequencies for /3 > —4.5 which is required from the observational constraints. The 
presented results are for non-extremal black holes, i.e., for bP^ < 1/8. Thus, the end of 
the sequences of black-hole solutions is reached when th ^>- which, for example, for 
M = 1 corresponds to charge P ^ 1.67. The results in the case when extremal black 
holes are present(i.e., for bP'^ > 1/8) are qualitatively the same. 

Before proceeding further, we will give some remarks on the computational abilities 
of the numerical method we use. As the numerical method is unstable in the region 
where the imaginary part of u is greater that the real part, we can compute only the 
fundamental [n = 0) QNM frequency when / = 0, because even for Schwarzschild black 
hole, the n = 1, / = mode is Un=i = 0.086 + 0.348z and the ratio of the imaginary 
part to the real part is beyond the computational abilities of the method we use. 

The results for the / = 2 fundamental modes as functions of the charge are shown 
on Fig. [2] for several values of the parameter (3 and for M = 1, and as we can see, 
both the real and the imaginary part of the modes increase as we increase the charge q 
The results are qualitatively the same for other values of / > 1. The positivity of 
the imaginary part according to our sign convention f ll7p means that the modes are 
damped and this is expected from the fact that the potential is positive for this values 
of the parameters. 

The situation changes when / = because in this case a negative minimum of the 
potential can exist when /3 < 0. On Fig. |3l the fundamental I = QNM frequencies 
are shown as a function of the charge for M = 1 and for several values of the parameter 
j3. Some of the calculated numerical values for the QNM frequencies are presented in 
Table [H The real and the imaginary part of the frequencies increase with the increase 
of the charge for (3 > 0. But for /3 < and for large values of the charge, both the real 
and the imaginary part of the frequencies start to decrease and eventually reach zero 
at some critical charge P^^-^^ q For larger values of the charge the real part is zero and 
the imaginary part is negative which means that the perturbations of the scalar field 
are growing with time and the corresponding black-hole solution is unstable. 

In the region where the real and the imaginary part of the frequency decreases, 
the imaginary part of the frequency is greater that the real one and the difference 



°For M = 1 and /3 = —4, the potential is positive for / > 1. 

^The critical charges where the real and the imaginary part reach zero coincide within 3% and the 
small discrepancy is due to the numerical instabilities described below. 
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Figure 2: The real ior (left panel) and the imaginary w/ (right panel) part of the / = 2 
fundamental (n = 0) frequencies as a function of the charge for several values of the 
parameter /3 and for M = 1.0. 



Table 1: 

Results 

given. 



The calculated numerical values of the / = fundamental QNM frequencies, 
for M = 1 and for several values of the parameter /3 and the charge P are 





(3 = 2 




/3 = 




(3 = -2 




/3 = -4 


p 


UjR Ui 




Ur Ui 




UjR UJi 




UJR Ui 


0.00 


0.1102 0.1049 




0.1102 0.1049 




0.1102 0.1049 




0.1102 0.1049 


0.20 


0.1111 0.1050 




0.1111 0.1050 




0.1110 0.1050 




0.1110 0.1050 


0.40 


0.1142 0.1056 




0.1135 0.1056 




0.1129 0.1055 




0.1122 0.1054 


0.60 


0.1211 0.1071 




0.1180 0.1067 




0.1148 0.1062 




0.1114 0.1058 


0.80 


0.1349 0.1106 




0.1252 0.1087 




0.1149 0.1068 




0.1035 0.1049 


1.00 


0.1613 0.1189 
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between them gets greater as we approach the critical charge P^J^ which means that 
in this region numerical instabilities are present. Indeed, the numerical results show 
that the computed values of the frequencies in this region depend on the values of the 
auxihary parameters, such as r^oo and r^^match, but the difference is up to 5%. Another 
justification of our results is the following. The critical value of the magnetic charge 
P^J^ can be calculated in two independent ways using the imaginary part of the modes 
in the stable and the unstable sectors because the method of calculation of the modes 
in each of them is different. Much to our pleasant surprise, the results for pI:J^^ obtained 
with the two independent methods coincide within 1%. 

The unstable QNM frequencies which are shown on Fig. [3l correspond only to the 
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Figure 3: The real ujr (left panel) and the imaginary ui (right panel) part of the / = 
fundamental QNM frequencies as a function of the magnetic charge P for several values 
of the parameter /3 and for M = 1 . 



first bound state of the potential well. With the increase of P, the potential well gets 
deeper, and the number of bound states (unstable modes) increases. We will label the 
different bound states with k, where fc = 0, 1, 2... . The fc-th unstable mode originates 
at -Pcj-it' where P^J^ < P^J^ < P^J^ < ..., and the radial part of the corresponding 
perturbation function has k zeros. The branches of the unstable modes of the first few 
bound states, are shown on Fig. H] for M = 1 and /3 = —4. 

It is clear that the first bound state (with A; = 0) is a continuation of the n = 
QNM frequencies into the unstable region on the uJi{P) diagram. We expect that the 
second bound state (with /c = 1) is a continuation of the n = 1 QNM frequencies and 
so on, but we can not check this numerically because as we said for / = only the 
fundamental n = mode can be calculated by the applied numerical method. 



5 Bifurcations and black hole non-uniqueness 

The presence of unstable modes (for I = 0) gives us a reason to expect that new 
spherically symmetric black-hole solutions with nontrivial scalar field exist and should 
bifurcate from the trivial solution. Some examples of such black-hole solutions have 
already been found in [20] . In the phase diagrams in [20], however, only one bifurcation 
point is present, and our numerical results show that it is just the point where the n = 
mode reaches zero lj. The scalar field of the non-trivial solution in [20] that bifurcates 
from the trivial solution is monotonous function of the radial coordinate, and when we 
take into account that the scalar field is zero at infinity, it is clear that it has no zeros. 
In the more general context of the present paper, it is natural to expect that solution 

(k) 

bifurcations occur for all of the critical charges P^j.^^ where branches of unstable modes 



^°In [5D], the charge P is kept fixed and the mass A/ is varied. With the decrease of M, the 
properties of the solutions change in the same way as in the case considered here when M is kept 
fixed and P is increased. 
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Figure 4: The imaginary part Ui as a function of the charge for the first four bound 
states {(jJr = in this case). 



originate. In other words, bifurcations should occur for the same values of the charge 
P for which the frequencies of the unstable modes become zero. The numbeio of 
the nontrivial solutions that bifurcate from the trivial solution should be equal to the 
number of the unstable modes - as we increase P, the number of bound states increases, 
and for each bound state that is added to the spectrum, a new bifurcation point occurs 
on the phase diagram at the same value of P. Actually, the static zero-modes (i.e., 
modes with w = and / = 0) are perturbative solutions of the field equations describing 
static, spherically symmetric black holes with non-trivial scalar field (i.e. hairy black 
holes) in the near vicinity of the bifurcation points. Consequently, we should expect 
that the scalar field of the non-trivial black-hole solutions originating at the bifurcation 
points should have the same number of zeros as a function of r as the corresponding 
unstable zero- mode. 

In order to check our hypothesis, we should solve the full system of field equations 
(fT2|) - flT^ . For this we will follow a modified form of the procedure from [20]. The 
system is coupled and nonlinear, and it has to be solved numerically. The domain of 
integration is from the radius of the horizon to infinity r G [th, oo) and the boundary 
conditions on both boundaries are 



firu) = 0, 



(35) 



lim m{r) 

r— j-oo 



M, 



(36) 
(37) 



lim 6{r) = lim (p{r) = 0, 

r—^oo r—^oo 

where M is the mass of the black hole in the Einstein frame. The following regulariza- 



^^We mean the number of nontrivial solutions up to the obvious discrete Z2 symmetry ip — > —tp of 
the field equations. 
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tion condition should also be fulfilled on the horizon 



f-^- Aa{^)A\^)[XdxLiX) - L(X)] 
ar ar 



(St 



r=rH 



The system (fT2|) - (fT4l) together with the boundary conditions is a BVP. The mass 
is an input parameter, and the radius of the horizon has to be determined additionally 
which means that the left boundary is a priori unknown. 

The BVP was solved numerically in [20] for the same class of scalar-tensor theories 
we consider in the present paper, and phases of the black-hole solutions were found 
when /3 < 0, i.e., in some regions of the parameter space, the number of black-hole 
solutions with the same mass and charge may be more than one, so an additional 
parameter must be introduced in order to label the different solutions. It is convenient 
to choose this parameter to be the scalar chargq^j defined by 



V 



lim r 

r— ^oo 



2^ 

dr 



(39) 
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Figure 5: The T>{P) dependence for branches of black-hole solutions for fixed value of 
the mass M = 1.0 and for /3 = —4. The trivial branch of solutions and the nontrivial 
branches characterized by scalar field which has up to three zeros are shown. The right 
panel is a magnification of the enclosed region of the 'D{P) diagram on the left panel. 



^^Another proper choice could be the value of scalar field on the horizon. Let us also make some 
clarifying remarks on the non-uniqueness of the black-hole solutions. It is twofold. For fixed values of 
the mass M and charge P the BVP described in this section has multiple solutions describing black 
holes with different scalar charges but also different radii of the event horizon. From mathematical 
point of view it would be unusual to interpret this situation as non-uniqueness since the different 
solutions have different domains [vh, oo). The presented formulation of the problem is more germane 
to the theory of black holes since the classification of solutions is usually based on the conserved charges 
such as the mass and electric or magnetic charge. However, another formulation of the problem is 
possible where the more usual for mathematics non-uniqueness of solution is observed. If we fix the 
charge P and the domain of integration [th, oo) we obtain black- hole solutions which have the same 
radius of the event horizon but different masses M. 



15 



M=l.0,fi=-4.0 



Trivial Solution 

Nontriv. Solution with zeros 

1 zero 

2 zeros 




M=1.0,/?=-4.0 -.^ 

Trivial Solution 

■ Nontriv. Solution with zeros 

1 zero 

- 2 zeros 




1.55 



1.60 



1.65 



1.70 



Figure 6: The radius of the horizon r// as a function of the charge P for some of the 
branches of black-hole solutions on Fig. |5l The right panel is a magnification of the 
enclosed region of the rniP) diagram on the left panel. 



Here, for the numerical solution of the above defined BVP, we apply two numerical 
methods in order to verify our results - the first method is the continuous analog of 
the Newtonian method that has been applied in {20] and the other one is the shooting 
method [26]. The solutions obtained with these two methods are in agreement. 

After performing a thorough search for new black-hole solutions, we indeed found 

(k) 

that the values of the charge P^^^^ where the unstable modes corresponding to the 
k-th bound state becomes zero and the values of the charge -Pbjfur corresponding to 
the k-th bifurcation point on the phase diagram T>{P) coincide. An example for this 
correspondence can be seen on Figs. |l]and[5]where the unstable modes and the branches 
of black-hole solutions are shown for values of the parameters M = 1 and /3 = —4. 



m 



m 



The charges P^^^^ and P^^f^^ are denoted on both figures, and for each value of k, they 
coincide within 0.2% error. On Fig. [Bjthe r^iP) dependence is shown for some of the 
branches of black-hole solutions presented on Fig. |5l 

As we said above, it is expected that the scalar field for the nontrivial solutions 
originating at -P^iifiir ^^^ ^ zeroes just as the perturbation function ^{r) for the unstable 
modes corresponding to the k-th bound state. The metric functions 6{r) and /(r) and 
the scalar field <f{r) of the nontrivial solutions for fixed values of the parameters M = 1, 
P = 1.66 and P = —4 are shown on Figs. [7] and [8] and as we can see indeed the k-th 
black-hole solution is characterized by a scalar field which has k zeros. An interesting 
observation is that the function /(r) for the nontrivial solutions with scalar field which 
has one or more zeros is not always monotonous and can have minima and maxima |lf|. 
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It can be proven that the function f(r) of the trivial solution is always monotonous. 
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Figure 7: The metric function 6{r) (left panel) and the scalar field y:>{r) (right panel) 
for nontrivial black-hole solutions for M = 1.0, P = 1.66, and /3 = —4. The solutions 
characterized by monotonous scalar field and by scalar field which has one and two 
zeros are shown. 
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Figure 8: The metric function /(r) for the same black-hole solutions as on Fig. [71 The 
right panel is a magnification of the enclosed region of the /(r) dependence on the left 
panel. 
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6 Stability analysis of the nontrivial solutions 

The next step in our investigations is to determine whether the obtained new nontrivial 
black-hole solutions are stable or not. Our approach is to study the QNMs associated 
with the pure radial perturbations of these solutions. As we will not consider general 
perturbations, only the instability of a black-hole solution can be proven rigorously 
within this approach. As we will see below, the nontrivial solutions for which the 
scalar field has zeros, are indeed unstable. 

The equations governing the radial perturbations of the nontrivial solutions were 
derived in [29] . It was shown in [29] that the radial perturbations of the metric and the 
electromagnetic field are determined by the perturbations of the scalar field. Hence, it 
is sufficient to consider only the perturbations of the scalar field 5ip which are governed 
by the equation 



VfV^''^''5ip - U{r)5ip = 0, 



(40) 



where 



U{r) 



l + 2rU'((/^)L(X(°))j \dMr 



n2 



+A\ip) [l6rdMr)aiip) + 49^«(^)] [x^'^^dx^L (x(°)) - L (x(°))] (41) 
-16a^(¥.)^^(¥.)[(xW)'4(o,L(xW) -XW9^(o,L(X(°)) +L(X(°))" 

and the covariant derivative V)( and the function X^°) are with respect to the un- 
perturbed background solution. After separating the variables and introducing the 
tortoise coordinate r* 



S^(r,t) = ^e'^\ 



dr^ 



dr 



we can obtain the following Schrodinger-like equation 

d'^^{r^) 



(42) 



dr'i 



+ UJ ipin) = f/eff(r*)V'(r* 



where the effective potential is given by 

f/es(r.) = /(r,)e-2^('-*){f/(r,) + 2A\^)L (X^) + 



(43) 



[l-/(r.)] . (44) 



The effective potential of the nontrivial black-hole solutions has negative minimum 
for all of the studied values of the parameters which means that unstable modes can 
exist. An example is given on Figs. |9] and [TO] where the effective potential Uesi'f') is 
plotted for the nontrivial solutions with scalar field which has no zeros and one or two 
zeros. 

We use the same numerical procedure for calculating the QNMs of the nontrivial 
solutions as the one described above. The numerical results show that for all of the 
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Figure 9: The potential Ucs of the nontrivial solutions characterized by monotonous 
scalar field (/^(r) for M = 1, /3 = —4 and for several values of the magnetic charge P. 
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Figure 10: The potential Ucs of the nontrivial solutions characterized by scalar field 
if{r) which has 1 zero (left panel) and 2 zeros (right panel) for M = 1, /3 = —4 and for 
several values of the magnetic charge P. 
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Figure 11: The imaginary part of the unstable modes (the real part is zero) corre- 
sponding to the first bound state of the trivial solutions and of the nontrivial solutions 
with scalar field ^p{r) which has up to three zeros (M = 1 and /3 = —4). 

studied values of the parameters the nontrivial solutions characterized by monotonous 
scalar field are stable against radial perturbations. The new nontrivial solutions char- 
acterized by scalar field which has one or more zeros are unstable. The lowest unstable 
modes as a function of the charge are shown on Fig. [TT] for the trivial solution and 
for the nontrivial solutions characterized by scalar field which has up to three zeros 
(the 'D{P) diagram of these black-hole solutions is given on Fig. [5]). The branches of 
unstable modes of the nontrivial solutions bifurcate from the trivial one at some critical 
charges which coincide with the charges -Pbifur where the nontrivial solutions bifurcate 
from the trivial one on the T>{P) diagram (Fig. [5]). The results also show that for all 
of the nontrivial solutions with k zeros, k bound states exist. 

The results about the stability of the nontrivial solutions remain valid for all of the 
studied values of the parameter. This means that all of the new nontrivial black hole 
branches characterized by scalar field fir) which has one or more zeros, are unstable, 
and the branches characterized by monotonous scalar field are stable against radial 
perturbations. The results in the previous section show that the trivial solution is 
stable against arbitrar}o (not only radial) perturbations before the first bifurcation 



>(o) 



(0) 



point -Pci-it- W^ expect that after the first bifurcation point -Pd-it' the solution possessing 
a scalar field without zeros is stable against arbitrary perturbations not only against 
the radial ones. 



7 Discussion 

In the present paper, we studied the scalar QNMs of the Einstein-Born-Infeld black- hole 
solution (the trivial solution) embedded in a certain class of scalar-tensor theories with 



^^The trivial solution is stable against general linear perturbations of the metric and electromagnetic 
field and this was demonstrated in 1501. 
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coupling function a{(p) = f3(p. The investigation of the spectrum of the QNMs shows 
that for /3 < 0, there are unstable modes which signal the presence of new nontrivial 
scalar-tensor, asymptotically flat, black- hole solutions with primary scalar hair that 
bifurcate from the trivial solution. These solutions were constructed numerically by 
solving the full system of static and spherically symmetric field equations. The number 
of nontrivial solutions that bifurcate from the trivial one is equal to the number of the 
bound states of the potential governing the scalar perturbations of the trivial solution. 
The nontrivial scalar-tensor black holes can be classified by the number of the zeros of 
the scalar field in the black hole exterior. We have shown that the scalar-tensor black 
holes possessing scalar field with one or more zeros are unstable. It seems that only 
the scalar-tensor black holes having scalar field without any zeros are stable solutions. 

The non-uniqueness of scalar-tensor black holes naturally raises the question of 
their classification. In other words, in the spirit of the no-hair conjecture, the question 
is whether some kind of classification theorem (conjecture) could be formulated. The 
most difficult question, however, is what kind of suitable further parameters associated 
with the black solutions should be specified in addition to the conserved asymptotic 
charges (i.e., the mass and the charge). The situation in some sense is similar to the 
situation in the higher dimensions (see, for example, [3I],[S2])- One possible solution 
of the problem is the non-conserved scalar charge V to be included in the classification 
theorem as one of the necessary additional parameters just as the dipole charge in the 
uniqueness theorem for the higher dimensional Einstein-Maxwell gravity [32]. Another 
obvious possibility is the inclusion of some "topological" characteristics of the scalar 
field (number of zeros, critical points) as additional parameters in the classification 
conjecture. Finally, we may adopt the view that the classical no-hair conjecture holds 
but only for the stable solutions. This point of view, however, seems to be not true. 
Our preliminary investigations show that, though for /3 < —4.5, probably there exists 
a domain in the parameter space where the trivial solution and the solution possessing 
scalar field without zeros are both linearly stable. In any case, the raised questions will 
need much deeper investigation. 

The results in the present paper also open a way for studying dynamical processes in 
the strong field regime within the framework of scalar-tensor theories. Especially, our 
results are the first steps to the study of strong field dynamical transitions from unstable 
to stable scalar-tensor black hole configurations. On the level of presented perturbative 
analysis, the characteristic time scales of the transitions between the unstable and 
stable black holes are of order r ~ l/w where u is the frequency of the corresponding 
lowest unstable mode. A deeper investigation of the dynamics of strong field transitions 
will require nonlinear analysis similar to that done in [33] for the scalar-tensor boson 
stars. 
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